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Abstract 

In this work the late-time evolution of Bianchi type VIII models 
is discussed. These cosmological models exhibit a chaotic behaviour 
towards the initial singularity and our investigations show that to- 
wards the future, far from the initial singularity, these models have 
a non-chaotic evolution, even in the case of vacuum and without in- 
flation. These space-time solutions turn out to exhibit a particular 
time asymmetry. On the other hand, investigations of the late-time 
behaviour of type VIII models by another author have the result 
that chaos continues for ever in the far future and that these solutions 
have a time symmetric behaviour: this result was obtained using the 
approximation methods of Belinski, Khalatnikov and Lifshitz (BKL) 
and we try to find out a possible reason explaining why the differ- 
ent approaches lead to distinct outcomes. It will be shown that, at a 
heuristic level, the BKL method gives a valid approximation of the 
late-time evolution of type VIII models, agreeing with the result of 
our investigations. 
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1 Introduction. 



The fundamental laws of physics are often said to be time symmetric: 
Newton's laws, Electromagnetism, Quantum Mechanics and Einstein's 
theory of gravity do not distinguish between past and future. However 
the phenomena taking place in our world display an arrow of time. 
It is of great interest to try to understand the physical origin of this 
time asymmetry. In this work we will focus on Einstein's theory of 
gravity and study how some solutions of the Einstein Field Equations 
(EFE) can distinguish between past and future. We will discuss the 
particular way in which the nature of the solutions at early times close 
to the initial singularity differs from the nature of the solution at late 
times, in the far future. 

The solutions that will be discussed are the Bianchi type VIII 
models, which are among the most general homogeneous solutions to 
the EFE. The type VIII model has been studied mainly at early 
times, where the possibility of chaotic behaviour was investigated. 
This solution was found by BKL [1],[2],[3] and Misner C. [4] to exhibit 
chaotic mixmaster behaviour as one approaches the initial singularity. 
As clock time i ^ the type VIII model evolves through an infinity of 
Kasner ( Bianchi type / ) stages, in a chaotic way, see [5] . The question 
arises to see whether this chaotic mixmaster behaviour continues for 
ever in the far future. 

P. Halpern [6] has used the methods of V. Belinski, I. M. Khalat- 
nikov and I. Lifshitz {BKL) [1] to study the late-time behaviour of 
type VIII models and the result is that chaos continues for ever in 
the far future of these models. Our investigation in [7], [8] shows that 
this can not be the case. The chaotic mixmaster solution close to the 
initial singularity evolves into a simpler non-chaotic solution. This 
result is consistent with the recent analysis of Ringstrom H. in [9] and 
other investigations supporting this result are [10], [11], [12] and [13]. 

It is important to understand why the application of the BKL 
method in [6] yields a completely different result. Either there could 
be some problem with the application to the late-time regime, or cither 
there could be something intrinsically wrong with the BKL method. 
BKL applied their method also to the study of the general inhomo- 
geneous solution at early times close to the initial singularity, so it is 
of relevance to determine what the nature of the problem actually is 
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and to see whether the BKL approximation could be used in the late- 
time regime as well. In order to achieve this we will compare in some 
detail the analysis of P. Halpern in [6] with the investigation of [7] , [8] . 

2 The time evolution of Bianchi type 
VIII models. 

The Bianchi type VIII universe is one of the four most general models 
of the spatially homogeneous universes and it is an eternally expanding 
solution. In recent work it is shown that type VIII universes become 
the most general of the spatially homogeneous models when space is 
compactified, see [16]. One of the reasons to believe that spatial sec- 
tions of cosmological models could be compact is that formulations of 
quantum cosmology require the volume of the spatial sections to be 
finite, in order for cosmological wave- functions to exist. So one could 
expect that it is most probable for a compact homogeneous universe 
to have type VIII symmetry. 

In [6] Halpern P. argues that the chaotic behaviour characteristic of 
the early time evolution of type VIII models continues for ever in the 
far future and that in this case there would be no essential difference 
between past and future for these models. However our investigation 
in [7], [8] shows that the behaviour in the past is very different from 
that in the future. 

Let us review how the Bianchi type VIII model evolves at late 

times, far from the initial singularity, according to [7], [8], [9], [10], 
[11], [12] and [13]. It turns out to be useful to discuss the problem 
in the Hamiltonian formalism [4]. In this formulation, using the so 
called metric approach [22] , where the metric components are the basic 
variables of the gravitational field, one can introduce group-invariant 
and time-independent frame vectors such that the line element for 
Bianchi class A models is given by, using an arbitrary time variable t' 

ds"^ = -N{t')dt'^ + gahWW^ 
where W"' are time-independent one-forms dual to the frame vec- 
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tors ea- We can introduce three time-dependent scale factors as 

gab = diag{a'^,h^,c^) 
which can be rewritten as 

gab = diag{e'^\e'P^e'^^) (1) 

with 

/3i = /30 - 2/3+ 

/?2 = /?o + /3+ + V^I3- 

Ps = + (3+ - V3p- 

The evolution of the type VIII model corresponds to the motion of 
a universe point in a triangular shaped potential in minisuperspace 
,P~), possessing an infinite open channel along the axis [22] . 
The motion of the universe point would then be analogous to motion 
of a ball in a triangular billiard, the difference being that the reflec- 
tion angle before the bounce with a wall is not equal to the angle 
after the bounce. If the universe point moves initially in a straight 
line ( corresponding to Kasner-like behaviour ) it will bounce off a 
potential wall and then move along another straight line ( a transi- 
tion to another Kasner-like behaviour ). The triangular shape of the 
potential allows the evolution to be chaotic. As the model evolves in 
the future, the Hamiltonian picture implies that the triangular poten- 
tial contracts, thus enabling the universe point to bounce of the walls 
in a chaotic sequence, ad infinitum, such that mixmastcr behaviour 
would continue for ever. Our analysis based on a combination of the 
Hamiltonian formalism and the orthonormal frame approach shows 
that this can not be the case. The universe point is forced to leave 
the triangular region of the potential and to escape along the infinite 
open channel along the axis, such that P~ and P'^ +00. 
This late-time evolution will be characterised by an infinity of non- 
chaotic oscillations between the two walls of the channel. The open 
channel becomes increasingly narrow as f3~^ ^ -|-oo and the universe 
point will exhibit increasingly rapid non-chaotic oscillations about the 
axisymmetric type VIII solution. The line element for vacuum type 
VIII models will tend in this way to the Bianchi type /// form of 
flat space-time. 
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ds"^ = -dt^ + [dx"^ + e^^dy^) + dz'^ (2) 

( the Bianchi type Vlh diagonal plane wave metric [22] with h = —1 
reduces to this line element for a particular choice of parameters ) such 
that the shear parameter S^, 
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( with cr^ being the shear scalar and H being the Hubble param- 
eter ), exhibits increasingly rapid non-chaotic oscillations. For non- 
vacuum perfect fluid models, with equation of state p = (7 — l)p 
and 7 > 2/3, the mixmaster behaviour will evolve as well to a simpler 
non-chaotic solution corresponding to the motion of the universe point 
along the infinite open channel of the potential. In this case the line 
element will tend to the Collins type /// [22] solution if2/3<7<l 
and to the above vacuum type /// line element if 1 < 7 < 2. In each 
of these cases the shear parameter exhibits increasingly rapid non- 
chaotic oscillations. Another remarkable property of these models is 
that the Weyl scalar [10], [12] is unbounded towards the far future: 
this might be of relevance in view of the relation between gravitional 
entropy and the Weyl curvature as conjectured by R. Penrose [19]. In 
this context, a further interesting study of the behaviour of the Weyl 
scalar for type VIII models and for other homogeneous solutions can 
be found in [20]. 



3 Is there a problem with the BKL 
approximation method ? 

Let us study the analysis of Halpcrn P. [6] and try to understand why 
this work leads to different results. 

First, in this work it is argued that in the Hamiltonian picture, 
the evolution of two neighbouring points in the type VIII potential is 
such that paths will diverge even if the points were close to each other 
initially. This is a feature of chaotic systems. Furthermore it is argued 
that this Hamiltonian system is time reversible and that therefore the 
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Bianchi type VIII model is chaotic both towards and away from the 
initial singularity. Our investigations show that the time reversibility 
can not be used to deduce the behaviour far from the singularity. The 
Bianchi type VIII system on the contrary exhibits time asymmetry, 
corresponding to the fact that as one approaches the singularity one 
has chaotic mixmaster behaviour while the far future is characterised 
by a simpler non-chaotic evolution. 

Second, another argument giving different results in the analysis 
of [6] is the following. In [6] essentially the same method was used as 
in the work of BKL [1], giving the result that chaos continues for ever 
in the far future. So docs this implies that there is a problem with 
the BKL method itself or has the method been applied incorrectly ? 
Or is the BKL approximation method not applicable in the late-time 
regime ? 

In [6] the field equations 

2(lna)" = (62 + c2)2-a^ (3) 
2(ln6)" = (a2 + c2)2-64 (4) 

2{lnc)" = {a'-b^f-c' (5) 

(Ina2)" + (ln62)" + (lnc2)" = {lrLa^y{lnb'^y+{lnby{lrLC^y+{lnc^y{lna^)' 

(6) 

(the prime denotes differentiation w.r.t. 77 defined by dt = {abc)dr]) 
were studied with the assumption that one of the scale factors is 
smaller then the other two, in order to study the evolution of a Kas- 
ner stage. The transition from one Kasner epoch to the other was 
parametrised by writing Kasner exponents as functions of a single 
parameter u 

u 

1 — u-\- 

1-u 
1 — u-\- 

— u 
1 — u + w' 

with u > 1 and such that the scale factors correspond to 
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Supposing that initially pi > > and p2 < 0, and if one includes 
only the dominant terms in (3)-(6), we get that, as t ^ +00, the field 
equations lead to another Kasner stage with p'3 > p'2 > and p\ < 
and so as the model evolves away from the singularity, the transition 
between Kasner epochs is given by the transformation 

+ 1 (7) 

Two of the scale factors will have exchanged their increasing and 
decreasing behaviour while the other scale factor continues to increase 
monotonically. Now, (7) means that the parameter u increases indef- 
initely, so it will never reach values u < 1, which implies that there is 
no reason for a new Kasner era to begin and the alternation between 
Kasner states seems to continue for ever. However, as the parameter 
u continues to increase, the Kasner parameters pi , p2 and ps approach 
the following values 

1 

u 

-1 

U 



and thus the exponents pi and p2 become close to each other. This 
means that the approximation done at the beginning of this calcula- 
tion, namely keeping the dominant scale factor a only in the right-hand 
side of (3)- (5), is not valid any more since a ^ b. In [6] the follow- 
ing notation is used to write the resulting field equations : a = Ina, 
P = \nb and 7 = Inc. When taking both a and b into account, (3)-(5) 
can be written as 

a +P =0 
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The equations (3)-(6) also yield a first integral which is given by, 
when terms containing c are neglected compared to those containing 
dominant contributions of a and 6, 

The subsequent analysis in [6] follows the one of [1]: the following 
change of time variable is made 



^oexp[^^(?7-?7o)] (8) 



with Q-o cind ryo being constants. If we define further % = a — /?, 
then the above field equations can be written as [6], [1] 

Xee + ^X? + ^sinh2x = (9) 

7^ = + + cosh2x + 1] (10) 

Now it was argued that as r/ — > +oo, we have ^ — > +oo because of 
equation (8), so that terms containing | can be neglected in (9) and 

(10 ). If we assume that a and h arc very close to each other, then we 
have sinh (2%) ~ 2% and these approximations lead to the following 
solution of (9)-(10) 

a, h = aoihhxp [±{^)sm (^ - Co)] (H) 

c = coexp[^(e2-^o')] (12) 

At some stage the increasing scale factor c will become compa- 
rable to a and b and will become the dominant term as ^ ^ +oo : 
this enables a new Kasner era to begin, with this time a,b. Fur- 
ther developments in [6] are done in order to show that this type of 
evolution is chaotic in the far future. 

Now, the solution (11)-(12) is valid if ^ — > -|-oo, which corresponds 
to T] ^ +00 because of equation (8). However as was explained in 
[7], a consequence of the inequality {abc\' > ( which can be deduced 
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directly from the field equations (3)-(5), the dot denotes diferentiation 
w.r.t. synchronous time t ) can be that 



V ^ (13) 

as t — > +CX). 

We can check this with the future asymptotic form of the line- 
element for vacuum Bianchi type VIII that we obtained, namely the 
Bianchi type /// form of flat space-time given by (2). For this line 
element we have 

abc 

This means that the relation between the ry-time and the syn- 
chronous time t is given by 
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dt -1 
abc t 
so that indeed 77 — > as t — -l-oo. 

This behaviour of BKL (rj) time given by equation (13) has also 
been observed for type Bianchi VIIq models in the work of Wainwright 
et al. [21]. 

This differs however with the derivation of BKL [1], ( see their 
equation 4.18 ) where the integration of dt = {abc)dr] yielded 

where ^0 is a constant, because this implies that ^ +00 as 
t — > +00. Such behaviour of ^ is also supposed to be true in the 
work of I. M. Khalatnikov and V. L. Pokrovsky [17]. This hypothesis 
is the source of the problem. Our result implies that ^ constant 
as t +00, which means that far from the initial singularity the 
solution (11)-(12) is no longer valid, since |-terms are not ignorable 
in equations (9)-(10). Indeed equation (9) with sinh(2x) 2% is a 
Bessel equation, a solution being a Bessel function of order zero 

xiO = MO 

The fact that ^ is assumed to continually increase is important in 
the analysis of [6] since this allows the scale factor c given by (12) to 
become larger then a and b, such that the intermediate axisymmetric 
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stage comes to an end in order for a new Kasner era to develop. As a 
result it was then explained in [6] that an " axisymmetric period" with 
two scale factors close to each other is finite in duration and that this 
ultimately has no effect on the solution. 

However if ^ — > constant as i — > +00, then xiO — ^oiO ~^ 
stant : the difference between a and /? will not change. In addition, 
the solution of equation (10) for 7(^) implies that 7(^) — constant 
as ^ ^ constant, such that the scale factor c can not become greater 
then a and b as was supposed in [6] . In other words the axisymmetric 
regime could never come to an end, which agrees with the result of 
our investigation and with that of other authors. 

Thus summarising one can say that provided one takes into account 
the correct interval for the BKL time variable r/, i.e. ] — 00, 0], then 
the BKL approximation method predicts that the vacuum Bianchi 
type VIII solution tends to an axisymmetric solution as t ^ +00. 
Thus as heuristic or approximative analysis the BKL method does 
not fail in the late-time regime. 

The rigorous analysis of Ringstrom H. in [18] of the Bianchi type 
IX model close to the singularity gave support to the validity of the 
BKL approximation method. In fact we have a new support to the 
heuristic validity of the BKL approximation method at late times, 
since it leads basically to a similar result as other rigorous methods 
[9], [10] and [12]. 

4 Chaos disappears in the far future. 

The fact that chaotic behaviour disappears at late times such that the 
type VIII models evolve into a simpler solution even when 7 > 2/3, 
thus without inflation, is remarkable in the sense that for most com- 
plex physical systems chaos is a feature of the whole evolution, for 
instance in the case of the three-body problem there seems to be no 
time direction in which chaos ceases to exist. In [15] it is explained 
that there is a frequently overlooked distinction between irreversible 
and chaotic behaviour of Hamiltonian systems. The latter does not 
in general appears to possess a direction of time, i.e. there is no es- 
sential distinction between past and future. By looking at a sequence 
of "snapshots" or configurations of the system at different instants 



10 



of time, one would not observe some asymmetry. The general solu- 
tions to the Einstein Field Equations that we discussed do not exhibit 
this distinction between chaos and time asymmetry, since the time 
evolution of the type VIII models distinguishes between past and fu- 
ture. The answer to the question "Is the type VIII system chaotic?" 
should then be: "It depends on the direction of time: towards the 
initial singularity the evolution is chaotic, while in the far future it is 
non-chaotic" . The type VIII system appears to loose memory of the 
chaotic initial state, and evolves in the far future to the axisymmetric 
type VIII solution. Thus the simple late-time asymptotic solution 
can be seen as arising from a wider set of ( chaotic ) initial conditions. 
This time asymmetry corresponds to the one defined by Halliwcll J. 
in [14], which is based on the distinction between the metric at early 
and late times. An important point to note is that the singularity of 
the type VIII solutions plays a crucial role in the dynamics of these 
models : some terms in the field equations become always negligeable 
(such that Kasner stages are always possible) close to the singularity, 
while far from the singularity they are not (see previous section) . This 
is related to the time asymmetry exhibited by the type VIII model. 

Finally, let tis point out that the asymmetry between past and fu- 
ture of type VIII models might be of relevance in the study of the 
late-time behaviour of cosmological solutions of (super)string theory 
as studied in [23] , where mixmaster behaviour was found to occur close 
to the initial singularity. 



5 Conclusion. 

In the work of BKL [1] and Halpcrn [6] the interval of variation of the 
BKL time variable r] is supposed to be ] — oo, +oc [. As a consequence 
the BKL approximation method predicts mixmaster behaviour to oc- 
cur in vacuum type VIII models at late times, far from the initial 
singularity. In this case the whole evolution of the type VIII model 
would be chaotic. If one takes into account that the correct interval 
of variation for r] is ] — oo,0[ then the BKL method predicts that 
at late times the type VIII model tends to the axisymmetric solu- 
tion, thus in agreement at a heuristic level with other rigorous meth- 
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ods. The BKL method is thus a vahd approximation at late times 
and it could be applied in order to obtain information about the na- 
ture of the general vacuum inhomogeneous solution in the late-time 
regime. The latc-timc evolution of type VIII models is non-chaotic, 
thus they distinguish between past and future. Remarkably the type 
VIII model looses memory of the chaotic initial state and evolves to a 
simpler non-chaotic future state. One could expect or conjecture that 
a similar time asymmetry will be found in the evolution of the general 
vacuum inhomogeneous solution. This remarkable type of dynamical 
feature of the Einstein Field Equations is not present in Newtonian 
or non-relativistic physics and is closely related to the occurence of 
singularities in these space-time solutions. 
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